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ON DEGENERATE PARTIAL DIFFERENTIAL EQUATIONS 

GUI-QIANG G. CHEN 

Abstract. Some of recent developments, including recent results, ideas, techniques, 
P^j and approaches, in the study of degenerate partial differential equations are surveyed 

and analyzed. Several examples of nonlinear degenerate, even mixed, partial differen- 
tial equations, are presented, which arise naturally in some longstanding, fundamental 
problems in fluid mechanics and differential geometry. The solution to these funda- 
mental problems greatly requires a deep understanding of nonlinear degenerate partial 
differential equations. Our emphasis is on exploring and/or developing unified math- 
ematical approaches, as well as new ideas and techniques. The potential approaches 
we have identified and/or developed through these examples include kinetic approaches, 
free boundary approaches, weak convergence approaches, and related nonlinear ideas and 
Qh ' techniques. We remark that most of the important problems for nonlinear degenerate 

•^T partial differential equations are truly challenging and still widely open, which require 

further new ideas, techniques, and approaches, and deserve our special attention and 
1 "T~i ' further efforts. 



1. Introduction 

> 

£T) ■ We survey and analyze some of recent developments, including recent results, ideas, 

techniques, and approaches, in the study of degenerate partial differential equations. We 
start with several important examples of degenerate/mixed linear degenerate equations 
and some of their interrelations. Then we present several examples of nonlinear degen- 

fundamental problems in fluid mechanics and differential geometry. These examples in- 
dicate that some of important nonlinear degenerate problems are ready to be tractable. 
Our emphasis is on exploring and/or developing unified mathematical approaches, as well 
as new ideas and techniques. The potential approaches we have identified and/or devel- 
oped through these examples include kinetic approaches, free boundary approaches, weak 
convergence approaches, and related nonlinear ideas and techniques. 

In fact, nonlinear degenerate, even mixed, partial differential equations arise also nat- 
urally in fundamental problems in many other areas such as elasticity, relativity, opti- 
mization, dynamical systems, complex analysis, and string theory. The solution to these 
fundamental problems in the areas greatly requires a deep understanding of nonlinear 
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degenerate partial differential equations. On the other hand, most of the important prob- 
lems for nonlinear degenerate partial differential equations are truly challenging and still 
widely open, which require further new ideas, techniques, and approaches, and deserve 
our special attention and further efforts. 

During the last half century, three different types of nonlinear partial differential equa- 
tions (elliptic, hyperbolic, parabolic) have been systematically studied separately, and 
great progress has been made through the efforts of several generations of mathemati- 
cians. With these achievements, it is the time to revisit and attack nonlinear degenerate, 
even mixed, partial differential equations with emphasis on exploring and developing uni- 
fied mathematical approaches, as well as new ideas and techniques. 

The organization of this paper is as follows. In Section 2, we present several important 
examples of linear degenerate, even mixed, equations and exhibit some of their interre- 
lations. In Section 3, we first reveal a natural connection between degenerate hyperbolic 
systems of conservation laws and the Euler-Poisson-Darboux equation through the en- 
tropy and Young measure, and then we discuss how this connection can be applied to 
solving hyperbolic systems of conservation laws with parabolic or hyperbolic degeneracy. 
In Section 4, we present a kinetic approach to handle a class of nonlinear degenerate 
parabolic-hyperbolic equations, the anisotropic degenerate diffusion-advection equation. 
In Section 5, we present two approaches through several examples to handle nonlinear 
mixed problems, especially nonlinear degenerate elliptic problems: free-boundary tech- 
niques and weak convergence methods. In Section 6, we discuss how the singular limits to 
nonlinear degenerate hyperbolic systems of conservation laws via weak convergence meth- 
ods can be achieved through an important limit problem: the vanishing viscosity limit 
problem for the Navier-Stokes equations to the isentropic Euler equations. 

2. Linear Degenerate Equations 

In this section, we present several important examples of linear degenerate, even mixed, 
equations and exhibit some of their interrelations. 

Three of the most basic types of partial differential equations are elliptic, parabolic, and 
hyperbolic. Consider the partial differential equations of second-order: 

d d 

y] <Hj (x)d XiX j u + ^bj (x)d Xi u + c(x)u = /(x), x = (x 1 ,x 2 ,--- ,Xd) € Q. C R d , (2.1) 

i,j=l j=l 

where ciy (x) , fy (x) , c(x) , and /(x) are bounded for x € Q. Equation (|2.1j) is called a 
uniformly elliptic equation in Q provided that there exists Ao > such that, for any 

£=(£i,-",&)eR d , 

d 
J^ <Hj(x)ti£j > Ao|£| 2 for x G ft, (2.2) 

that is, the d x d matrix (ojj(x)) is positive definite. 

Two of the basic types of time-dependent partial differential equations of second-order 
are hyperbolic equations: 

d tt u-L x u = f(t,x.), t>0, (2.3) 
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and parabolic equations: 

d t u-L x u = f(t,x), t>0, (2.4) 

where L x is a second-order elliptic operator: 

d d 

L*u= } ] Ojjjt, x)d XiXl u + } y bj(t, x)d x .u + c(t, x)u, (2.5) 

i,j=l j=l 

for which Oij(t, x),bj(t, x),c(£,x), and /(t,x) are locally bounded for (£, x) € M + + := 
[0, oo) x M rf . Equations (|2.3p and (|2.4p are called uniformly hyperbolic or parabolic equa- 
tions, respectively, in a domain in M + + under consideration, provided that there exists 
A > such that, for any £ £ R d , 

d 

^Mi,x)^>A |£| 2 (2.6) 

on the domain, that is, the dxd matrix (ay(i,x)) is positive definite. In particular, when 
the coefficient functions a^, 6j, c, and the nonhomogeneous term / are time-invariant, and 
the solution u is also time- invariant, then equations (|2.3|) and (|2.4|) are coincide with 
equation (|2.ip . 

Their representatives are Laplace's equation: 

A x u = 0, (2.7) 

the wave equation: 

d a u - A x u = 0, (2.8) 

and the heat equation: 

d t u - A x n = 0, (2.9) 

respectively, where A x = ^j=i ®xjxj is the Laplace operator. 

Similarly, a system of partial differential equations of first-order in one-dimension: 

d t U + A(t,x)d x U = 0, Pel" (2.10) 

is called a strictly hyperbolic system, provided that the n x n matrix A(t, x) has n real, 
distinct eigenvalues: 

Xl(t,x) < A 2 (t, x) < ■ ■ ■ < \ n (t,x) 

in the domain under consideration. The simplest example is 

A(t,x)=(_° 1 - 1 ). (2.11) 

Then system ()2.10p - (|2.1ip with U = (u,v) T is equivalent to the two wave equations: 

dttu - d xx u = 0, d u v - d xx v = 0. 

However, many important partial differential equations are degenerate or mixed. That 
is, for the linear case, the matrix (ajj(t, x)) or (ajj(x)) is not positive definite or even 
indefinite, and the eigenvalues of the matrix A(t, x) are not distinct. 
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2.1. Degenerate Equations and Mixed Equations. Two prototypes of linear degen- 
erate equations are the Euler-Poisson-Darboux equation: 

(x — y)d xy u + ad x u + f3d y u = 0, (2-12) 

or 

x(d xx u- dyyu) +ad x u + (3d y u = 0, (2-13) 

and the Beltrami equation: 

x(d xx u + dyyu) + ad x u + /3d y u = 0. (2-14) 

Three prototypes of mixed hyperbolic-elliptic equations are the Lavrentyev-Betsadze 
equation: 

d xx u + sign(x)dyyU = (2-15) 

which exhibits the jump from the hyperbolic phase x < to the elliptic phase x > 0, the 
Tricomi equation: 

d xx u + xd yy u = (2.16) 

which exhibits hyperbolic degeneracy at x = (i.e., two eigenvalues coincide at x = 0, 
but the corresponding characteristic curves are not tangential to the line x = 0), and the 
Keldysh equation: 

xd xx u + dyyU = (2-17) 

which exhibits parabolic degeneracy at x = (i.e., two eigenvalues coincide at x = 0, but 
the corresponding characteristic curves are tangential to the line x = 0). 

The mixed parabolic-hyperbolic equations include the linear degenerate diffusion-advection 
equation: 

d t u + b- Vu = V- (AVu), xGlR^, t e M+ := (0,oo), (2.18) 

where u : M + x IR — >• M is unknown, V = (<9 Xl , . . . , 9 Xd ) is the gradient operator with 
respect to x = (xi, ...,x d ) € R d , and b : R + x M d ->■ R d and A : M + x R d -> R dxd 
are given functions such that A = (a,ij(t, x)) is a d x d nonnegative, symmetric matrix. 
When the diffusion matrix function A degenerates, the advection term b • Vu dominates; 
otherwise, the parabolic diffusion term V • (A Vu) dominates. 

2.2. Interrelations between the Linear Equations. The above linear equations are 
not actually independent, but are closely interrelated. 

2.2.1. The Wave Equation (|2.8p and the Euler-Poisson-Darboux Equation (|2.13j) . Seek 
spherically symmetric solutions of the wave equation (I2.8P : 



v(t,r) = n(t,x), r = |x|. 

Then v(t, r) is governed by the Euler-Poisson-Darboux equation (J2.13P with a = d — 1 and 
(3 = 0: 

duv — d rr v d r v = 0. 

r 

As is well-known, the Euler-Poisson-Darboux equation plays an important role in the 

spherical mean method for the wave equation to obtain the explicit representation of the 

solution in R d , d > 2 (cf. Evans [69]). 



d TT U + OyyU + — d T u = 0, 
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2.2.2. The Tricomi Equation (|2.16p . the Beltrami Equation (|2.14p . and £/ie Euler-F 'oisson- 
Darboux Equation (|2.13p . Under the coordinate transformations: 

0,y) — > (r,y) = (|(-x)2,y), when x < 0, 

(x,y) — )• (T,y) = (fx2,y), when x > 0, 

the Tricomi equation (|2.16p is transformed into the Beltrami equation (|2.14j) with a = ^ 
and /3 = when x = r > 0: 

1 

3r" 
and the Euler-Poisson-Darboux equation (|2.13p with a = i and /3 = when x = r < 0: 

d TT u — d vv u H <9 T u = 0. 

yy 3r 

These show that a solution to one of them implies the solution of the other correspondingly, 
which are equivalent correspondingly. 

Linear degenerate partial differential equations have been relatively better understood 
since 1950. The study of nonlinear partial differential equations has been focused mainly 
on the equations of single type during the last half century. The three different types of 
nonlinear partial differential equations have been systematically studied separately, and 
one of the main focuses has been on the tools, techniques, and approaches to understand 
different properties and features of solutions of the equations with these three different 
types. Great progress has been made through the efforts of several generations of mathe- 
maticians. 

As we will see through several examples below, nonlinear degenerate, even mixed, par- 
tial differential equations naturally arise in some fundamental problems in fluid mechanics 
and differential geometry. The examples include nonlinear degenerate hyperbolic systems 
of conservation laws, the nonlinear degenerate diffusion-advection equation and the Euler 
equations for compressible flow in fluid mechanics, and the Gauss-Codazzi system for iso- 
metric embedding in differential geometry. Such degenerate, or mixed, equations naturally 
arise also in fundamental problems in many other areas including elasticity, relativity, op- 
timization, dynamical systems, complex analysis, and string theory. The solution to these 
fundamental problems in the areas greatly requires a deep understanding of nonlinear 
degenerate, even mixed, partial differential equations. 

3. Nonlinear Degenerate Hyperbolic Systems of Conservation Laws 

Nonlinear hyperbolic systems of conservation laws in one-dimension take the following 
form: 

d t U + d x F(U) = 0, U£R n , (t,x) G M+. (3.1) 

For any C 1 solutions, (13. If) is equivalent to 

d t U + VF(U)d x U = 0, U e R n , (t, x) e R 2 + . 

Such a system is hyperbolic if the n x n matrix VF(U) has n real eigenvalues \j(U) and 
linearly independent eigenvectors Tj{U), 1 < j < n. Denote 

V := {U : A,(t0 = Xj{U), i?j,l<i,j< n} (3.2) 
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as the degenerate set. If the set T> is empty, then this system is strictly hyperbolic and, 
otherwise, degenerate hyperbolic. Such a set allows a degree of interaction, or nonlinear 
resonance, among different characteristic modes, which is missing in the strictly hyperbolic 
case but causes additional analytic difficulties. A point £/* G T> is hyperbolic degenerate 
if VF(U#) is diagonalizable and, otherwise, is parabolic degenerate. 

Degenerate hyperbolic systems of conservation laws have arisen from many important 
fields such as continuum mechanics including the vacuum problem, multiphase flows in 
porous media, MHD, and elasticity. On the other hand, degenerate hyperbolicity of sys- 
tems is generic in some sense. For example, for three-dimensional hyperbolic systems 
of conservation laws, Lax |102| indicated that systems with 2(modA) equations must be 
degenerate hyperbolic. The result is also true when the systems have ±2, ±3, ±4(mod8) 
equations (see [75] ) . Then the plane wave solutions of such systems are governed by the 
corresponding one-dimensional hyperbolic systems with T> ^ 0. 

Since F is a nonlinear function, solutions of the Cauchy problem for (|3.ip (even starting 
from smooth initial data) generally develop singularities in a finite time, and then the 
solutions become discontinuous functions. This situation reflects in part the physical 
phenomenon of breaking of waves and development of shock waves. For this reason, 
attention is focused on solutions in the space of discontinuous functions, where one can 
not directly use the classical analytic techniques that predominate in the theory of partial 
differential equations of other types. 

To overcome this difficulty, a natural idea is to construct approximate solutions U e (t,x) 
to (|3.ip by using shock capturing methods and then to study the convergence and con- 
sistency of the approximate solutions to (13. ID . The key issue is whether the approximate 
solutions converge in an appropriate topology and the limit function is consistent with 
(13. ip . Solving this issue involves two aspects: one is to construct good approximate solu- 
tions, and the other is to make suitable compactness analysis in an appropriate topology. 

3.1. Connection with the Euler-Poisson-Darboux Equation: Entropy and Young 
Measure. The connection between degenerate hyperbolic systems of conservation laws 
and the Euler-Poisson-Darboux equation is through the entropy. 

A pair of functions (rj, q) : R n — > M 2 is called an entropy-entropy flux pair (entropy pair, 
in short) if they satisfy the following linear hyperbolic system: 

Vq(U) = Vr)(U)VF(U). (3.3) 

Then the function rj{U) is called an entropy. Clearly, any C 1 solution satisfies 

d t r](U) + V7](U)VF(U)d x U = 0, (3.4) 

or 

dtv(U) + d x q(U) = 0. (3.5) 

For a BV solution that is not C , the second term in fj3.4|) has no meaning in the classical 
sense because of the multiplication of a Radon measure with a discontinuous function. If 
r](U) is an entropy, then the left side of (|3.4p becomes the left side of (j3.5[) that makes 
sense even for L°° solutions in the sense of distributions. An L°° function U(t, x) is called 
an entropy solution (cf. Lax |101| ) if 

d t 7](U) + d x q(U) < (3.6) 

in the sense of distributions for any convex entropy pair, that is, the Hessian V 2 r/(C7) > 0. 
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Assume that system (|3.ip is endowed with globally defined Riemann invariants w(JJ) = 
(wi, . . . , w n )(U), 1 < i < n, which satisfy 

Vwi(U) ■ VF(U) = \i(U)Vwi(U). (3.7) 

The necessary and sufficient condition for the existence of the Riemann invariants Wi(U), 1 < 
i < n, for strictly hyperbolic systems is the well-known Frobenius condition: 

li{rj,r k } = 0, for any j, k / i, 

where L, denotes the left eigenvector corresponding to Aj and {•, ■} is the Poisson bracket of 
vector fields in the [/-space (cf. [U [99]). For n = 2, the Frobenius condition always holds. 
For any smooth solution U(t,x), the corresponding Riemann invariants itfj, 1 < i < n, 
satisfy the transport equations 

d t Wi + Xi(w)d x Wi = 0, 

which indicate that Wi is invariant along the i-th characteristic field. 

Taking the inner product of (j3.7|) with the right eigenvectors r, of VF produces the 
characteristic form 

(AiVr/ - V?) ■ r t = 0, 
that is, 

K0 Wi i] = d w .q. (3.8) 

For n = 2, the linear hyperbolic system (|3.7p is equivalent to (|3.8p . and system (|3.8|) 
can be reduced to the following linear second-order hyperbolic equation: 

d wi \2 n d W2 Xi a , . 

o WlW2 V + T ^Ow 2 V ~ 7 ^o Wl r] = 0. (3.9) 

A2 — M A2 — Ai 



For the case T> = 0, equation (|3.9p is, in general, regular and hyperbolic for which either 
the Goursat problem or the Cauchy problem is well posed in the coordinates of Riemann 
invariants. The entropy space is infinite-dimensional and is represented by two families 
of functions of one variable. However, for the case T> ^ 0, the situation is much more 
complicated because of the singularity of the functions A "j|_ A 2 and ^ry- on the set w(£>) C 
M 2 . The typical form of such equations is 

a(^, Wl -w 2 ) P(^,wi-w 2 ) 
d U!1 w2 r l^ d wi r]^ 2 - d W2 r] = 0. (3.10) 

W\ — W2 Wi — W2 

This equation may have extra singularities of a and (3 both at the origin and on the line 
w\ = W2- The questions are whether there exist nontrivial regular solutions of the singular 
equation (|3.10p and, if so, how large the set of smooth regular solutions is. 

The connection between the compactness problem of approximate solutions and the 
entropy determined by the Euler-Poisson-Darboux equation (|3.9p is the Young measure 
via the compensated compactness ideas first developed by Tartar [1471 1148] and Murat 
|119| 1120] and a related observation presented by Ball [7] . 

The Young measure is a useful tool for studying the limiting behavior of measurable 
function sequences. For an arbitrary sequence of measurable maps 

[7 e :lRi^M n , lllHIroo <C< oo, 
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that converges in the weak-star topology of L°° to a function U, 

w* - lim U e = U, 

there exist a subsequence (still labeled) U e and a family of Young measures 
v t ,x € Prob.(M n ), supp v t ,x C {A : |A| < M}, (t,x)eM+, 

such that, for any continuous function g, 



w* - lim g(U e (t,x)) = / g(\)dv t JX) := (u ttX ,g) 

for almost all points (t,x) G R^_. In particular, U e (t,x) converges strongly to U(t,x) if 
and only if the Young measure Vt,x is equal to a Dirac mass concentrated at U(t, x) for 
a.e. (t, x). 

In many cases, one can estimate (cf. |61|, I30|) that the approximate solutions U e (t,x) 
generated by the shock capturing methods for (J3.3J) satisfy 

• II^IIl 00 < c < co. 

• For C 2 entropy pairs (r]i,qi),i = 1,2, determined by (|3.8p and (|3.9h . 

d t ^(C/ e ) + ^%(C/ e ) compact in fl^J. (3.11) 

Then, for any C 2 entropy pairs (r]i,qi),i = 1,2, determined by f|3.8j) and (|3.9p . the Young 
measure ut jX is forced to satisfy the Tartar-Murat commutator relation: 

(vt,x,mq2- qim) = (vt,x,m){vt,x,q2) - (vt,x,qi){vt,x,m) a.e. (t,x)eR+. (3.12) 

This relation is derived from ()3.11[) . the Young measure representation theorem for the 
measurable function sequence, and a basic continuity theorem for the 2x2 determinant 
in the weak topology (cf. [1471 [30]). This indicates that proving the compactness of 
the approximate solutions U e (t,x) is equivalent to solving the functional equation (|3.12p 
for the Young measure vt,x f° r a H possible C 2 entropy pairs (r]i,qi),i = 1,2, determined 
by (|3.8p and ()3.9p . If one clarifies the structure of the Young measure f t)X (A), one can 
understand the limiting behavior of the approximate solutions. For example, if one can 
prove that almost all Young measures Uf jX , (t, x) G M. 2 ,, are Dirac masses, then one can 
conclude the strong convergence of the approximate solutions U e (t,x) almost everywhere. 
On the other hand, the structure of the Young measures is determined by the C 2 solutions 
of the Euler-Poisson-Darboux equation (|3.9p via (|3.12p . One of the principal difficulties for 
the reduction is the general lack of enough classes of entropy functions that can be verified 
to satisfy certain weak compactness conditions. This is due to possible singularities of 
entropy functions near the regions with degenerate hyperbolicity. The larger the set of C 2 
solutions to the entropy equation (|3.9p of the Euler-Poisson-Darboux type is, the easier 
solving the functional equation (I3.12p for the Young measure is. 

3.2. Hyperbolic Conservation Laws with Parabolic Degeneracy. One of the pro- 
totypes for such systems is the system of isentropic Euler equations with the form 

(d t p + d x m = 0, 

\d t m + d x (f+p( P )) = o, [ ' ; 

where p, m, and p are the density, mass, and pressure, respectively. For p > 0, v = m/p 
represents the velocity of the fluid. The physical region for (I3.13P is {(/?, m) : \m\ < C p} 
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2 

for some C > 0, in which the term — in the flux function is only at most Lipschitz 
continuous near the vacuum. For p > 0, v = — represents the velocity of the fluid. 
The pressure p is a function of the density through the internal energy e(p): 

p(p) = p 2 e'(p) forp>0. (3.14) 

In particular, for a polytropic perfect gas, 



p(r) 



p(p) = K P \ e(j>)= ^dr = -?—pi-\ (3.15) 

Jo r 7-1 

where 7 > 1 is the adiabatic exponent and, by the scaling, the constant k in the pressure- 
density relation may be chosen as k = 4 without loss of generality. 

For (|3.13p . strict hyperbolicity and genuine nonlinearity away from the vacuum require 
that 

p'(p) > 0, 2p'(p) + pp"(p) > for p > 0. (3.16) 

Near the vacuum, for some 7 > 1, 

P(p) 



P< 



«i > whenp^O. (3.17) 



One of the fundamental features of this system is that strict hyperbolicity fails when 
p — > 0. That is, T> = {p = 0}, which is the vacuum state, and the degeneracy is parabolic. 
The eigenvalues of (|3.13p are 

Ai = c(p), A 2 = \-c(p), 

P P 

and the corresponding Riemann invariants of f|3. 13|) are 

in in 

w\ = k(p), w 2 = \-k(p), 

P P 

where 

c(p) := Vp'(p) 

is the sound speed, and 

Up) : = J" 2^1 dr . 
Jo r 
Then the entropy equation (|3.9p in the Riemann coordinates is the following Euler-Poisson- 
Darboux equation: 

d WlW2 r]-\ (dwiV - 9 W2 ri) = 0, (3.18) 

^1 — W2 



where 



withp = -A;- 1 (Hi^2; 
For the 7-law gas, 



k(p)k"(p) 

3-7 



a{w\ — W2) = X :- 



2(7-1) 
is a constant, the simplest case. 
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In terms of the variables (p,v), r\ is determined by the following second-order linear 
wave equation: 

d ppV -k'(p) 2 d vv ri = 0. (3.19) 

An entropy r\ is called a weak entropy if 77(0, v) = 0. For example, the mechanical energy 
pair: 

lm 2 m 3 , p(p), 

V* = o — + P e (/>> Q* = 5-5 + m{e{p) + J, 

2 p 2p z p 

is a convex weak entropy pair. 

By definition, the weak entropy kernel is the solution x{P-> v > s ) °f the singular Cauchy 

problem 

d PP X-k'(p) 2 d vv x = 0, 

X (0,v,a) = 0, (3.20) 

d p x(0,v,s) = 5 V=S 

in the sense of distributions, where s plays the role of a parameter and 5 V=S denotes the 
Dirac measure at v = s. Then the family of weak entropy functions is described by 

v(P,v)= / x(p,v,s)ip(s)ds, (3.21) 

Jm. 

where ip(v) is an arbitrary function. By construction, rj(0,v) = 0, rj p (0,v) = ip(v). One 
can prove that, for < p < C, |^| < C, 

\Vr](p,m)\ < C v , \V 2 r](p,m)\ < C^-q^p.m), 

for any weak entropy r/, with C v independent of (p,m). 

Since this system is a prototype in mathematical fluid dynamics, the mathematical 
study of this system has an extensive history dating back to the work of Riemann 



where a special Cauchy problem, so-called Riemann problem, was solved. Zhang-Guo 
|157| established an existence theorem of global solutions to this system for a class of 
initial value functions by using the characteristic method. Nishida [121j obtained the first 
large data existence theorem with locally finite total-variation for the case 7 = 1 using 
Glimm's scheme [80J. Large-data theorems have also been obtained for general 7 > 1 in 
the case where the initial value functions with locally finite total-variation are restricted to 
prevent the development of cavities (e.g. [60 | I108|. 1122] ). The difficult point in bounding 
the total-variation norm at low densities is that the coupling between the characteristic 
fields increases as the density decreases. This difficulty is a reflection of the fact that the 
strict hyperbolicity fails at the vacuum: D/0. 

The first global existence for (J3.13P with large initial data in L°° was established in 
DiPerna [6b for the case 7 = — ^-, N > 5 odd, by the vanishing viscosity method. The 
existence problem for usual gases with general values 7 G (1, |] was solved in Chen [29] and 
Ding-Chen-Luo [59] (also see [30]). The case 7 > 3 was treated by Lions- Perthame-Tadmor 
|107j . Lions-Perthame-Souganidis |106j dealt with the interval (|,3) and simplified the 
proof for the whole interval. A compactness framework has been established even for the 
general pressure law in Chen-LeFloch [44] by using only weak entropy pairs. 

More precisely, assume the pressure function p = p(p) G C 4 (0,oo) satisfies condition 
(I3.16P (i.e., strict hyperbolicity and genuine nonlinearity) away from the vacuum and, near 
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the vacuum, p(p) is only asymptotic to the 7-law pressure (as real gases): there exist a 
sequence of exponents 

1 < 7 : = 7l < 72 < . . . < lN < (37 - l)/2 < 7iV+ i (3.22) 

and a sufficiently smooth function P = P(p) such that 

JV 

p(p) = X>iP 7j '+P 7iV+ip (P)> ( 3 - 23 ) 

i=i 

P(p) and p 3 P'"(p) are bounded as p -> 0, (3.24) 

for some coefficients «,- £ M with «i > 0. The solutions under consideration will remain 
in a bounded subset of {p > 0} so that the behavior of p(p) for large p is irrelevant. This 
means that the pressure law p(p) has the same singularity as ^7=1 K j P lj near the vacuum. 
Observe that p(0) = p'(0) = 0, but, for k > 7 i, the higher derivative p^ k \p) is unbounded 
near the vacuum with different orders of singularity. 

Theorem 3.1 (Chen-LeFloch |44|). Consider system ()3.13|) to£/i general pressure law 
satisfying ()3.16|) and (|3.22p ()3.24j) . Assume that a sequence of functions (p £ ,m £ ) satisfies 
that 

(i) There exists C > independent of e such that 

0<p e {t,x)<C, \m £ (t,x)\<Cp £ (t,x) for a.e. (t,x)\ (3.25) 

(ii) For any weak entropy pair (rj,q) of (|3.13|) . (|3.16|) . and (|3.22j) - ([3.24j) . 

d t r](p £ ,m £ ) + d x q(p e ,m e ) is compact in H^l(R+). (3.26) 

Then the sequence (p £ ,m £ ) is compact in L i 1 oc (lR^_). Moreover, there exists a global en- 
tropy solution (p(t,x),m(t,x)) of the Cauchy problem (|3.13p . (|3.16p . and (|3.22|) - (|3.24|) . 
satisfying 

< p(t, x) < C, \m(t, x) I < Cp(t, x) 

for some C depending only on Co and 7 , and 

d t r)(p, m) + d x q(p, m) < 

in the sense of distributions for any convex weak entropy pair (r],q). Furthermore, the 
bounded solution operator (p,m)(t, •) = St(po,mo)(-) is compact in L 1 for t > 0. 

As discussed in Section 3.1, if the Young measure satisfying f|3. 12|) reduces to a Dirac 
mass for a.e. (t,x), then the sequence (p £ ,m e ) is compact in the strong topology and 
converges subsequentially toward an entropy solution. For the Euler equations, to obtain 
that the Young measure V(t, x ) ls a Dirac mass in the (p, m)-plane, it suffices to prove that 
the measure in the (p, u)-plane, still denoted by V(t,x)i i s either a single point or a subset 
of the vacuum line 



i m £ 



{(p,v) : p = 0,\v\< sup — - }. 

£>0 P £ UL 

The main difficulty is that only weak entropy pairs can be used, because of the vacuum 
problem. 

In the proof of [29], (59J [61] (also cf. [30]), the heart of the matter is to construct special 
functions ip in (|3.2ip in order to exploit the form of the set of constraints (|3.12p . These 
test-functions are suitable approximations of high-order derivatives of the Dirac measure. 
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It is used that (|3.12p represents an imbalance of regularity: the operator on the left is more 
regular than the one on the right due to cancellation. DiPerna [61J considered the case 
that A := 2 , ~^A > 3 is an integer so that the weak entropies are polynomial functions of 
the Riemann invariants. The new idea of applying the technique of fractional derivatives 
was first introduced in Chen [29J and Ding-Chen-Luo [59J to deal with real values of A. 

A new technique for equation (|3.12p was introduced by Lions-Perthame-Tadmor |107| 
for 7 G [3, oo) and by Lions-Perthame-Souganidis [106J for 7 € (1,3). Motivated by a 
kinetic formulation, they made the observation that the use of the test-functions ip could 
in fact be bypassed, and (|3.12p be directly expressed with the entropy kernel x*- Namely, 
(|3.12p holds for all si and S2 by replacing rjj := x*i s j) an d qj ■= o~*(sj) for j = 1, 2. Here 
cr* is the entropy-flux kernel defined as 

o-*(p, v,s) = (v + 6(s- v)) x*(p, v, s). 

In |106| . the commutator relations were differentiated in s, using the technique of fractional 
derivatives as introduced in |29(, 159] by performing the operator 9 A+1 , so that singularities 
arise by differentiation of x* ■ This approach relies again on the lack of balance in regularity 
of the two sides of (I3.12J) and on the observation that {v(t )X )iX*{ s )) is smoother than the 
kernel x*( s ) itself, due to the averaging by the Young measure. 

However, many of the previous arguments do not carry over to the general pressure 
law. The main issue is to construct all of the weak entropy pairs of (I3.13P . The proof in 
Chen-LeFloch [44] has been based on the existence and regularity of the entropy kernel 
that generates the family of weak entropies via solving a new Euler-Poisson-Darboux 
equation, which is highly singular when the density of the fluid vanishes. New properties 
of cancellation of singularities in combinations of the entropy kernel and the associated 
entropy-flux kernel have been found. In particular, a new multiple-term expansion has 
been introduced based on the Bessel functions with suitable exponents, and the optimal 
assumption required on the singular behavior on the pressure law at the vacuum has been 
identified in order to valid the multiple-term expansion and to establish the existence 
theory. The results cover, as a special example, the density-pressure law p(p) = K\ p 11 + 
K2 p 72 where 71,72 € (1,3) and «i,«2 > are arbitrary constants. The proof of the 
reduction theorem for Young measure has also further simplified the proof known for the 
polytropic perfect gas. 

Then this compactness framework has been successfully applied to proving the con- 
vergence of the Lax-Friedrichs scheme, the Godunov scheme, and the artificial viscosity 
method for the isentropic Euler equations with the general pressure law. 

3.3. Hyperbolic Conservation Laws with Hyperbolic Degeneracy. One of the 

prototypes of hyperbolic conservation laws with hyperbolic degeneracy is the gradient 
quadratic flux system, which is umbilic degeneracy, the most singular case: 

d t U + d x {VuC{U)) = 0, U = (u,v) T eR 2 , (3.27) 

where 

C(U) = -au 3 + bu 2 v + uv 2 , (3.28) 

and a and b are two real parameters. 

Such systems are generic in the following sense. For any smooth nonlinear flux function, 
take its Taylor expansion about the isolated umibilic point. The first three terms including 
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the quadratic terms determine the local behavior of the hyperbolic singularity near the 
umbilic point. The hyperbolic degeneracy enables us to eliminate the linear term by a 
coordinate transformation to obtain the system with a homogeneous quadratic polynomial 
flux. Such a polynomial flux contains some inessential scaling parameters. There is a 
nonsingular linear coordinate transformation to transform the above system into (I3.27h - 
(|3.28p . first studied by Marchesin, Isaacson, Plohr, and Temple, and in a more satisfactory 
form by Schaeffer-Shearer [135] . From the viewpoint of group theory, such a reduction from 
six to two parameters is natural: For the six dimensional space of quadratic mappings acted 
by the four dimensional group GL{2, R), one expects the generic orbit to have codimension 
two. 

The Riemann solutions for such systems were discussed by Isaacson, Marchesin, Paes- 
Leme, Plohr, Schaeffer, Shearer, Temple, and others (cf. [88j ESI H36| I140J ). Two new 
types of shock waves, the overcompressive shock and the undercompressive shock, were 
discovered, which are quite different from the gas dynamical shock. The overcompressive 
shock can be easily understood using the Lax entropy condition [101j . It is known that 
there is a traveling wave solution connecting the left to right state of the undercompres- 
sive shock for the artificial viscosity system. Stability of such traveling waves for the 
overcompressive shock and the undercompressive shock has been studied (cf. |109|, 1110] ). 

The next issues are whether the compactness of the corresponding approximate solutions 
is affected by the viscosity matrix as the viscosity parameter goes to zero, to understand the 
sensitivity of the undercompressive shock with respect to the viscosity matrix, and whether 
the corresponding global existence of entropy solutions can be obtained as a corollary from 
this effort. The global existence of entropy solutions to the Cauchy problem for a special 
case of such quadratic flux systems was solved by Kan [93] via the viscosity method. A 
different proof was given independently to the same problem by Lu [111] . 

In Chen-Kan [41) . an L°° compactness framework has been established for sequences 
of approximate solutions to general hyperbolic systems with umbilic degeneracy specially 
including (J3.27[) ( f3.28[) . Under this framework, any approximate solution sequence, which 
is apriori uniformly bounded in L°° and produces the correct entropy dissipation, leads to 
the compactness of the corresponding Riemann invariant sequence. This means that the 
viscosity matrix does not affect the compactness of the corresponding uniformly bounded 
Riemann invariant sequence. Again, one of the principal difficulties associated with such 
systems is the general lack of enough classes of entropy functions that can be verified 
to satisfy certain weak compactness conditions. This is due to possible singularities of 
entropy functions near the regions of degenerate hyperbolicity. The analysis leading to 
the compactness involves two steps: 

In the first step, we have constructed regular entropy functions governed by a highly 
singular entropy equation, the Euler-Poisson-Darboux type equation, 



! \ mo / <-v " v ', 



d WlW2 v + 2 d WlV + W2 d W2 r) = 0. (3.29) 

There are two main difficulties. The first is that the coefficients of the entropy equation are 
multiple-valued functions near the umbilic points in the Riemann invariant coordinates. 
This difficulty has been overcome by a detailed analysis of the singularities of the Rie- 
mann function of the entropy equation. This analysis involves a study of a corresponding 
Euler-Poisson-Darboux equation and requires very complicated estimates and calculations. 
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An appropriate choice of Goursat data leads to the cancellation of singularities and the 
achievement of regular entropies in the Riemann invariant coordinates. The second diffi- 
culty is that the nonlinear correspondence between the physical state coordinates and the 
Riemann invariant coordinates is, in general, irregular. A regular entropy function in the 
Riemann invariant coordinates is usually no longer regular in the physical coordinates. 
We have overcome this by a detailed analysis of the correspondence between these two 
coordinates. 

In the second step, we have studied the structure of the Young measure associated with 
the approximate sequence and have proved that the support of the Young measures lies 
in finite isolated points or separate lines in the Riemann invariant coordinates. This has 
been achieved by a delicate use of Serre's technique |137j and regular entropy functions, 
constructed in the first step, in the Tartar-Murat commutator equation for Young measures 
associated with the approximate solution sequence. 

This compactness framework has been successfully applied to proving the convergence 
of the Lax-Friedrichs scheme, the Godunov scheme, and the viscosity method for the 
quadratic flux systems. Some corresponding existence theorems of global entropy solutions 
for such systems have been established. The compactness has been achieved by reducing 
the support of the corresponding Young measures to a Dirac mass in the physical space. 



4. Nonlinear Degenerate Parabolic-Hyperbolic Equations 

One of the most important examples of nonlinear degenerate parabolic-hyperbolic equa- 
tions is the nonlinear degenerate diffusion-advection equation of second-order with the 
form: 

8feit + V-f(«) = V-(A(tt)V«), x€R d , teR+, (4.1) 

where f : R — > M. d satisfies f'(-) G Lg > c (R;R d ), and the d x d matrix A(u) is symmetric, 
nonnegative, and locally bounded, so that A(u) = (crik{u))(<Jik(u)) T with <Jik{u) G L^ C (R). 
Equation ()4.ip and its variants model anisotropic degenerate diffusion-advection motions 
of ideal fluids and arise in a wide variety of important applications, including two-phase 
flows in porous media and sedimentation-consolidation processes (cf. [191 123 E7J 1150] ). 
Because of its importance in applications, there is a large literature for the design and 
analysis of various numerical methods to calculate the solutions of equation (|4.ip and its 
variants (e.g. [271 ED E51 E7J EL] ) for which a well-posedness theory is in great demand. 

One of the prototypes is the porous medium equation: 

d t u = A x (u m ), m>l, (4.2) 

which describes the fluid flow through porous media (cf. |150| ). Equation (J4.2J) is degen- 
erate on the level set {u = 0}; away from this set, i.e., on {u > 0}, the equation is strictly 
parabolic. Although equation (J4.2J) is of parabolic nature, the solutions exhibit certain 
hyperbolic feature, which results from the degeneracy. One striking family of solutions is 
Barenblatt's solutions found in [8]: 



u(t, x; cl,t) 



1 -.2 



(t + r) 



k(m — 1) |x 
a~ 



2 , i 

771 — 1 



2md 



(t + r) 



d 
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where k = -, — _,,,, , and a/0 and r > are any constants. The dynamic boundary of 
the support of u(t, x; a, r) is 



2md N k 

k(m — 1) 



|x|=4/— — a(t + r)d, t>0. 



This shows that the support of Barenblatt's solutions propagates with a finite speed! 
The simplest example for the isotropic case (i.e., A (it) is diagonal) with both phases is 

u 2 
d t u + d x {—) = d xx [u]+, (4.3) 

where [u]+ = max{it,0} (cf. |35|). Equation fj4.3|) is hyperbolic when u < and parabolic 
when u > 0, and the level set {u = 0} is a free boundary that is an interface separating 
the hyperbolic phase from the parabolic phase. For any constant states u such that 
u + < < u~ and y* £ R, there exists a unique nonincreasing profile <ft = 4>{y) such that 

lim (p(y) = u ± , lim <f>'{y) = 0, 

y— ¥±oo y— >±oo 

4>{y) = u + for all y > y*, 

G C 2 (-oo,y*), 0(y*-) = O, lim <j)'(y) = -oo, 



dim\- 



i + - 

= —u^u 
y=y*- 2 



so that u(t,x) = <f>(x — st) is a discontinuous solution connecting u + from n _ , where 
s = - \ u . Although u(t,x) is discontinuous, [u(t,x)]+ is a continuous function even 
across the interface {u = 0}. See Fig. 1. 



Figure 1. A discontinuous profile connecting the hyperbolic phase from 
the parabolic phase 

The well-posedness issue for the Cauchy problem is relatively well understood if one re- 
moves the diffusion term V- (A(u)Vtt), thereby obtaining a scalar hyperbolic conservation 
law; see Lax [lQOj . Oleinik |123j . Volpert [151j . Kruzhkov [98] . and Lions- Perthame-Tadmor 
[107] . and Perthame [126] (also cf. [54 ], I143J ) . It is equally well understood for the diffusion- 
dominated case, especially when the set {u : rank(A(u)) < d} consists of only isolated 
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points with certain order of degeneracy; see Brezis-Crandall [17], Caffarelli- Friedman |21j . 
Daskalopoulos-Hamilton [55], DiBenedetto [56], Gilding [79], Vazquez |150j . and the ref- 
erences cited therein. For the isotropic diffusion, Oij(u) = 0, i ^ j, some stability results 
for entropy solutions were obtained for BV solutions by Volpert-Hudjaev |152] in 1969. 
Only in 1999, Carrillo [23] extended this result to L°° solutions (also see Eymard-Gallouet- 
Herbin- Michel [72] . Karlsen-Risebro [91j for further extensions), and Chen-DiBenedetto 
|35j handled the case of unbounded entropy solutions which may grow when |x| is large. 

A unified approach, the kinetic approach, to deal with both parabolic and hyperbolic 
phases for the general anisotropic case with solutions in L 1 has been developed in Chen- 
Perthame [35]. This approach is motivated by the macroscopic closure procedure of the 
Boltzmann equation in kinetic theory, the hydrodynamic limit of large particle systems 
in statistical mechanics, and early works on kinetic schemes to calculate shock waves and 
theoretical kinetic formulation for the pure hyperbolic case; see [HI [Ml [26j [961 H07[ 11251 
126, I127}ll44j and the references cited therein. In particular, a notion of kinetic solutions 
and a corresponding kinetic formulation have been extended. More precisely, consider the 
kinetic function, quasi-Maxwellian, x onl 2 : 

{+1 for < v < u, 
-1 for u < v < 0, (4.4) 

otherwise. 

If ue L°°(R + ;L l (R d )), then x(v,u) £ L°°(R + ; L l (R d+l )). 

Definition. A function u(t,x) G L°°(lR + ; L 1 (M d )) is called a kinetic solution if u(t, x) 
satisfies the following: 
(i) The kinetic equation: 

d t x{v;u) +?{v) -X7x(v,u) - V • (A(v)Vx(v; u)) = d v {m + n)(t,x,v) (4.5) 

holds in the sense of distributions with the initial data x( v 'i u )\t=o = x( v j u o), for some 
nonnegative measures m(t,x;v) and n(t,x;v), where n(t,x,v) is defined by 

d d „ 

(n(i,x,.),V(-)> :=E(E^tM*> x ))) ^(Rf 1 ), (4-6) 

fc=l i=l 



for any tp E C$°(R) with ip > and 0f k {u) := /" a ik (v)y/ip(v) dv; 
(ii) There exists \x G L°°(R) with < /j,(v) — )• as \v\ — > oo such that 

j I (m + n)(t,x;v)dtdx< fj,(v); (4.7) 

JO JM. d 

(iii) For any two nonnegative functions ip\, ip2 G Cq°(M), 

^Mu(t,x))J2d x ^t k 2 (u(t,x)) = ^Td Xi ^(u(t,x)) a.e. (t,x). (4.8) 

i=l i=l 

Then we have 

• Well-posedness in L : Under this notion, the space L 1 is both a well-posed space 
for kinetic solutions and a well-defined space for the kinetic equation in (i). That 
is, the advantage of this notion is that the kinetic equation is well defined even 
when the macroscopic fluxes and diffusion matrices are not locally integrable so 
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that L 1 is a natural space on which the kinetic solutions are posed. This notion also 
covers the so-called renormalized solutions used in the context of scalar hyperbolic 
conservation laws by Benilan-Carrillo-Wittbold |llj . Based on this notion, a new 
approach has been developed in [46] to establish the contraction property of kinetic 
solutions in L . This leads to a well-posedness theory — existence, stability, and 
uniqueness — for the Cauchy problem for kinetic solutions in L . 

• Consistency of the kinetic equation with the original macroscopic equa- 
tion: When the kinetic solution u is in L°°, for any n G C 2 with rj"(u) > 0, 
multiplying rf(v) both sides of the kinetic equation in (i) and then integrating in 
wet yield 

dtv( u ) + V x • (q(n) - A(u)S7 x r)(u)) = — rj"(v)(m + n)(t, x; v)dv < 0. 

Jr 

In particular, taking r](u) = ±w yields that u is a weak solution to the macroscopic 
equation. The uniqueness result implies that any kinetic solution in L°° must be 
an entropy solution. On the other hand, any entropy solution is actually a kinetic 
solution. Therefore, the two notions are equivalent for solutions in L°°, although 
the notion of kinetic solutions is more general. 

• For the isotropic case, condition (iii) automatically holds, which is actually a chain 
rule. In fact, the extension from the isotropic to anisotropic case is not a purely 
technical issue, and the fundamental and natural chain-rule type property (iii) 
does not appear in the isotropic case and turns out to be the corner-stone for the 
uniqueness in the anisotropic case. Moreover, condition (ii) implies that m + n has 
no support at u = oo. 

Based on this notion, the corresponding kinetic formulation, and the uniqueness proof in 
the pure hyperbolic case in |126j . we have developed a new effective approach to establish 
the contraction property of kinetic solutions in L , This leads to a well-posedness theory 
for the Cauchy problem of (I3.13P with initial data: 



u\ t =o = n G L°°(R d ) (4.9) 



for kinetic solutions only in L 1 . 



Theorem 4.1 (Chen-Perthame [46]). (i) For any kinetic solution u G L°°(1R + ; L^iW 1 )) 
with initial data uq{x), we have 

\\u(t,-)-u \\ L i (R d)^-0 as t-s>0; 



(ii) Ifu,v G L 00 (M_|_;L 1 (M d )) are kinetic solutions to (|4.ip and (|4.9p with initial data no( x ) 
andvo(x), respectively, then 

\\u(t, •) - v(t, 0Hli(r<i) < IN - v Q \\ LHR dy, (4.10) 



(iii) For initial data Uq G L 1 (M d ), there exists a unique kinetic solution u G C(R+; L l ( 
for the Cauchy problem A3. 131) and (14. 9p . If Uq G L°° n L x (IR d ), then the kinetic solution 
is the unique entropy solution and \\u(t, ■)||i oo (E <i ) — II^oIIl 
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Furthermore, assume that the flux function f (u) £ C 1 and the diffusion matrix A(n) £ C 
satisfy the nonlinearity-diffusivity condition: The set 

{v : t + f» • y = 0, yA(«) y T = 0}cl (4.11) 

has zero Lebesgue measure, for any r £ M, y = (j/i, . . . , j/d) lozf/i |y| = 1. Lei it £ L°°(IR + + ) 
be the unique entropy solution to (J3.13P and (|4,9p with periodic initial data uq £ L°° /or 
periodTp = Uf =1 [0,Pi\, i.e., uo( x + -Pj e i) = ^o( x ) o,.e., where {ei}f =1 is the standard basis 
inR d . Then 

1 f 

\\u(t,-)-—— u (x)dxL lfT , -»0 ast^oo, (4.12) 

where |Tp| is i/ie volume of the period Tp. 

The nonlinearity-diffusivity condition implies that there is no interval of t> in which 
the flux function f(y) is affine and the diffusion matrix A(t>) is degenerate. Unlike the 
pure hyperbolic case, equation (|4.ip is no longer self-similar invariant, and the diffusion 
term in the equation significantly affects the behavior of solutions; thus the argument in 
Chen-Frid [40] based on the self-similar scaling for the pure hyperbolic case could not be 
directly applied. The argument for proving Theorem 14.11 is based on the kinetic approach 
developed in [46], involves a time-scaling and a monotonicity-in-time of entropy solutions, 
and employs the advantages of the kinetic equation (|4.5|) , in order to recognize the role of 
nonlinearity-diffusivity of equation (14. 1|) (see [47] ) . 

Based on the very construction of the kinetic approach, the results can easily be trans- 
lated in terms of the old Kruzkov entropies by integrating in v. In the case of uniqueness 
for the general case, this was performed in Bendahmane-Karlsen [10]. 

Follow-up results based on the Chen-Perthame approach in [46] include L 1 -error esti- 
mates and continuous dependence of solutions in the convection function and the diffusion 
matrix in Chen-Karlsen [42]; and more general degenerate diffusion-advection-reaction 
equations in Chen-Karlsen [43]. Other recent developments include the related notion 
of dissipative solutions in Perthame-Souganidis |128j . regularity results of solutions in 
Tadmor-Tao [146j , as well as different types of diffusion terms in Andreianov-Bendahmane- 
Karlsen [5]. 

5. Nonlinear Equations of Mixed Hyperbolic-Elliptic Type 

Unlike the linear case, very often, nonlinear equations can not be separated as two 
degenerate equations, but are truly mixed. Nonlinear partial differential equations of 
mixed hyperbolic-elliptic type arise naturally in many fundamental problems. In this 
section, we present two approaches through several examples to handle nonlinear mixed 
problems, especially nonlinear degenerate elliptic problems: Free-boundary techniques and 
weak convergence methods. 

5.1. Weak Convergence Methods. We first present two problems: transonic flow pass 
obstacles and isometric embedding, for which weak convergence methods, especially meth- 
ods of compensated compactness, play an important role. 
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5.1.1. Transonic Flow Pass Obstacles in Gas Dynamics. By scaling, the Euler equations 
for compressible, isentropic, irrotational fluids take the form: 

(d x (pu) + d y (pv) = 0, 
\d x v-d y u = 0, 

combined with the Bernoulli relation: 

7-1/. 2 , 2n\^FT 



p={l-^ r (u' + v')y-\ 7 >1, (5.2) 

for the pressure-density relation: p = p(p) = p 7 /7- This provides two equations for the two 
unknowns (u,v). Furthermore, we note that, if p is constant (which is the incompressible 
case), the two equations in (|5.ip become the Cauchy-Riemann equations for which any 
boundary value problem can be posed for the elliptic partial differential equations. 
By the second equation in (|5.ip . we introduce the velocity potential (p: 

(u,v)=V<p. (5.3) 

Then our conservation law of mass becomes a nonlinear partial differential equation of 
second-order for ip: 

d x (pd x ip) + dy{pdyip) = 0, (5.4) 

which is combined with the Bernoulli relation: 

1-7-^lV^I 2 )^. (5.5) 

Introduce the sound speed c as 

c = p'(p) = 1 — q , with the fluid speed q = \Vp\ = \/u 2 + v 2 , (5.6) 

so that, at the sonic value when q = c, we have q = q* with 



1* := \l^i- ( 5 - 7 ) 



Then equation (|5.4p is elliptic if 

q < q* 
and hyperbolic when 

q > q*- 
There is an upper bound placed on q from the Bernoulli relation: 



q < q C av ■= y^TY' ( 5 - 8 ^ 

where q cav is the cavitation speed for which p = 0. 

On the other hand, equation (|5.4p corresponds to the Euler-Lagrange equation for the 
functional 



f G(\Vp\)dxdy, (5.9) 



where 

G(q)= J (J-— --s) '<!*. (5.1U) 



7-1 \ ^i 
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Since (g* — q)G"{q) > 0, the direct method of calculus of variations (e.g. Evans |69|) 
provides the existence of weak solutions if it is known apriori that the flow is subsonic 
(q < q*) so that G is convex and the problem is elliptic. For example, this includes the 
fundamental problem of subsonic flow around a profile as formulated in Bers [13] . 

A profile V is prescribed by a smooth curve, except for a trailing edge with an opening 
£7r at zt, < e < 1. 




Figure 2. Profile C of the obstacle V: e = 




Figure 3. Profile C of the obstacle V: e = 1 




Figure 4. Profile C of the obstacle V: < e < 1 

If e = 0, the profile has a tangent at the trailing edge (see Fig. [2|). The tangent to V 
satisfies a uniform Holder condition with respect to the arc length. The velocity w = (u, v) 
must attain a given subsonic limit at infinity. We enforce the Kutta-Joukowski condition: 

q — >■ as (x, y) — > zt if £ = 1, 

g = O(l) as (x, y) ->■ zt if < £ < 1, 

and define problem Pi(woo) with a prescribed constant velocity w^ at infinity. 
For a smooth profile, e = 1 (see Fig. [3]), define the circulation 



P 



(u,v) ■ t ds, 



(5.11) 

■ r) 



where r is the unit tangent to V. In this case, we can also consider problem P2( w c 
where the circulation is prescribed, instead of the Kutta-Joukowski condition. 
In both problems, the slip boundary condition on the profile is required: 

(u,v) ■ n = on V (boundary condition) (5-12) 

where n is the exterior unit normal on V. 
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The first existence theorem for P^Woo) was given by Frankl-Keldysh [74J for suffi- 
ciently small speed at infinity. For a general gas, the first complete existence theorem for 
■fb(woo,r) was given by Shiffman |141j . This was followed by a complete existence and 
uniqueness theorem by Bers [12] for Pi(w 0O ), a stronger uniqueness result of Finn-Gilbarg 
[73j . and a higher dimensional result of Dong [64 j. The basic result is as follows: 



For a given constant velocity at infinity, there exists a number q depending upon the 
profile V and the adiabatic exponent 7 > 1 such that the problem Pi(w 0O ) has a unique 
solution for < q^ := |woo| < q. The maximum q m of q takes on all values between and 
Q*, Qm — > as qoc — > 0, and q m — > q* as qoo — > q. A similar result holds for ^(woo, T). 

The main tool for the results is to know apriori that, if g^ < q (i.e., the speed at 
infinity is not only subsonic) but sufficiently subsonic, then q < q* in the whole flow 
domain. Subsonic flow at infinity itself does not guarantee that the flow remains subsonic, 
since the profile produces flow orthogonal to the original flow direction. Shiftman's proof 
did use the direct method of the calculus of variations, while Bers's relied on both elliptic 
methods and the theory of pseudo-analytic functions. The existence of a critical point for 
the variational problem would be a natural goal for the case when q^ is not restricted 
to be less than q, since it would provide a direct proof of our boundary value problem. 
However, no such proof has been given. 

More recent investigations based on weak convergence methods start in the 1980's. 
DiPerna [63] suggested that the Murat-Tartar method of compensated compactness be 
amenable to flows which exhibit both elliptic and hyperbolic regimes, and investigated 
an asymptotic approximation to our system called the steady transonic small disturbance 
equation (TSD). He proved that, if a list of assumptions were satisfied (which then guaran- 
teed the applicability of the method of compensated compactness), then a weak solution 
exists to the steady TSD equation. However, no one has ever been able to show that 
DiPerna's list is indeed satisfied. 

In |115j (also see |117j ). Morawetz layed out a program for proving the existence of 
the steady transonic flow problem about a bump profile in the upper half plane (which 
is equivalent to a symmetric profile in the whole plane). In particular, Morawetz showed 
that, if the key hypotheses of the method of compensated compactness could be satisfied, 
now known as a "compactness framework" (see Chen [21]), then indeed there would exist 
a weak solution to the problem of flow over a bump which is exhibited by subsonic and 
supersonic regimes, i.e., transonic flow. 

The "compactness framework" for our system can be stated as follows: A sequence of 
functions w e (x,y) = (u e ,v £ )(x,y) defined on an open set Oct 2 satisfies the following 
set of conditions: 

(A.l) q e (x,y) = |w e (x,y)| < g* a.e. in Q for some positive constant q* < q ca v', 

(A. 2) d x Qi±('w e )+dyQ2±('w £ ) are confined in a compact set in H^ c (Q) for entropy pairs 
(Qi±> Q2±)i and (Qi±(w £ ), <32±(w £ )) are confined to a bounded set uniformly in L^ C (Q), 
where (Qi,Q2) is an entropy pair, that is, d x Q\(w £ ) + d y Q2(~w £ ) = along smooth 
solutions of our system). 
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In case (A.l) and (A. 2) hold, then the Young measure v X y determined by the uniformly 
bounded sequence of functions w £ (x,y) is constrained by the following commutator rela- 
tion: 

\Vx,y> Q1+Q2- — Q1-Q2+) = {Vx,y, Ql+)( ly x,y, <?2-) ~ \ v x,y, Ql-)\^x,y, Q2+) ■ (5.13) 

The main point for the compensated compactness framework is to prove that u xy is 
a Dirac measure by using entropy pairs, which implies the compactness of the sequence 
w £ (x,y) = (u £ ,v £ )(x,y) in Lj oc (fi). In this context, both DiPerna [63] and Morawetz [115j 
needed to presume the existence of an approximating sequence parameterized by e to their 
problems satisfying (A.l) and (A. 2) so that they could exploit the commutator relation 
and obtain the strong convergence in L} (Q) to a weak solution of their problems. 

As it turns out, there is one problem where (A.l) and (A. 2) hold trivially, i.e., the sonic 
limit of subsonic flows. In that case, we return to the result by Bers [12] and Shiffman 
|141j . which says that, if the speed at infinity g^ is less than some q, there is a smooth 
unique solution to problems Pi(woc) and ^2(^00, T) and ask what happens as Qoo /*■ q. In 
this case, the flow develops sonic points and the governing equations become degenerate 
elliptic. Thus, if we set e = q — c/oo and examine a sequence of exact smooth solutions to 
our system, we see trivially that (A.l) is satisfied since \q £ \ < q*, and (A. 2) is satisfied 
since d x Q±(w £ ) + d y Q±(w £ ) = along our solution sequence. The effort is in finding 
entropy pairs which can guarantee that the Young measure v xy reduces to a Dirac mass. 
Ironically, the original conservation equations of momentum in fact provide two sets of 
entropy pairs, while the irrotationality and mass conservation equations provide another 
two sets. This observation has been explored in detail in Chen-Dafermos-Slemrod-Wang 



What then about the fully transonic problem of flow past an obstacle or bump where 
qoo > <P- In Chen-Slemrod-Wang [48J, we have provided some of the ingredients for 
satisfying (A.l) and (A. 2). More precisely, we have introduced the usual flow angle 6 = 
tan _1 (-) and written the irrotationality and mass conservation equation as an artificially 
viscous problem: 



d x v £ - d y u £ = eA6 £ , 

d x (p £ u £ ) + d y (p £ v £ ) = eV ■ (a(p £ )Vp 



v - (514) 



where c(p) is suitably chosen, and appropriate boundary conditions are imposed for this 
regularized "viscous" problem. The crucial point is that a uniformly L°° bound in q £ 
has been obtained when 1 < 7 < 3 which uniformly prevents cavitation. However, in 
this formulation, a uniform bound in the flow angle £ and lower bound of q £ away from 
zero (stagnation) in any fixed region disjoint from the profile are still assumed apriori to 
guarantee the (q, #)-version of (A.l). On the other hand, (A. 2) is easily obtained from 
the viscous formulation by using a special entropy pair of Osher-Hafez- Whitlow [124]. In 
fact, this entropy pair is very important: It guarantees that the inviscid limit of the above 
viscous system satisfies a physically meaningful "entropy" condition (Theorem 2 of [124] ). 
With (A.l) and (A. 2) satisfied, Morawetz's theory [115] then applies to yield the strong 
convergence in Lj (fi) of our approximating sequence. It would be interesting to establish 
a uniform bound in the flow angle 9 £ and lower bound of q £ away from zero (stagnation) 
in any fixed region disjoint from the profile under some natural conditions on the profile. 
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5.1.2. Isometric Embeddings in Differential Geometry. In differential geometry, a long- 
standing, fundamental problem is to characterize intrinsic metrics on a two-dimensional 
surface Ai 2 which can be realized as embeddings into R 3 . Let C R 2 be an open set and 
gij,i,j = 1,2, be a given matric on M 2 parameterized on Q. Then the first fundamental 
form for Ai 2 is 

I = gn(dx) 2 + 2g 12 dxdy + g 2 2(dy) 2 ■ (5.15) 

Isometric Embedding Problem. Seek an injective map r : — > R 3 such that 
dr ■ dr = I, i.e., 

d x r-d x r = g u , d x r ■ d y r = g 12 , d y r ■ d y r = g 22 , (5.16) 

so that (d x r, d y r) in R 3 are linearly independent. 

The equations above are three nonlinear partial differential equations for the three 
components of r. Recall that the second fundamental form 77 for Ai 2 defined on $7 is 

II = -dn ■ dr = h n (dx) 2 + 2h u dxdy + h 22 (dy) 2 , (5.17) 

and hij is the orthogonality of the unit normal n of the surface r(S7) C R 3 to its tangent 
plane. The Christoffel symbols are 

rf-^idjgu+digji-dwij), 

which depend on the first derivatives of gij, where the summation convention is used, (g ) 
denotes the inverse of (gij), and (di,d 2 ) := (d x ,d y ). 

The fundamental theorem of surface theory indicates that there exists a surface in R 3 
whose first and second fundamental forms are I and II if the coefficients (gij) and (hij) of 
the two given quadratic forms I and //, I being positive definite, satisfy the Gauss-Codazzi 
system. This theorem also holds even for discontinuous coefficients hij (cf. Mardare |114j ). 

Given (gij), the coefficients (hij) of II are determined by the Gauss-Codazzi system, 
which consists of the Codazzi equations: 



d x M - d y L = rg } L - 2vf 2 ] M + vf^N, 
d x N - d y M = -T^L + 2rf 2 } M - r^JV, 



(5.18) 



and the Gauss equation, i.e., Monge- Ampere constraint: 

LN-M 2 = K, (5.19) 

where 

L = -^, M=— p=, N = ^7Fji Isl = det(g i3 ) = g n g 22 - g\ 2 , (5.20) 

K(x, y) is the Gauss curvature that is determined by the relation: 

K(x,y) = ?f^, (5.21) 

\9\ 

and 

ttijkl - gim \p k l {j - C>jL ik + 1 tj 1 nk - 1 ik 1 nj ) 

is the Riemann curvature tensor depending on g^ and its first and second derivatives. 
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Therefore, given a positive definite (gij), the Gauss-Codazzi system consists of the three 
equations for the three unknowns (L, M, N) determining the second fundamental form. 
Note that, while (g%j) is positive definite, -R1212 may have any sign and hence K may have 
any sign. 

From the viewpoint of geometry, the Gauss equation is a Monge- Ampere equation and 
the Codazzi equations are as integrability relations. On the other hand, we are interested in 
a fluid mechanic formulation for the isometric embedding problem so that the problem may 
be solved via the approaches that have shown to be useful in fluid mechanics, especially for 
the nonlinear partial differential equations of mixed hyperbolic-elliptic type. To achieve 
this, the way to reformulate the problem via solvability of the Codazzi equations under 
the Gauss equation has been adopted in Chen-Slemrod-Wang [39]. 

Set 

L = pv 2 + p, M = —puv, N = pu 2 
choose the "pressure" p as for the Chaplygin type gas: p = 



+ P, 

— -, and set q 2 
Then the Codazzi equations become the balance equations of momentum: 



u 2 + v 2 . 



rf 2 W+p) 



2T 12 puv 



r { S(pu 2 + P ), 



d x (puv) + d y (pv +p) = - 

d x { P u 2 +p) + d y (puv) = -r^j(pv 2 + P ) - 2T^puv - r^ipu 2 + P ), 

and the Gauss equation becomes the Bernoulli law: 



1 



P 



Vq T TK' 



(5.22) 



(5.23) 



Set the "sound" speed: c = \/p'(p), i.e., c 2 



7 1 



q 2 + K. Then 



When K > 0, 
When K < 0, 
When K = 0, 



7 2 < c 2 and the "flow" is subsonic; 
7 2 > c 2 and the "flow" is supersonic; 
a 2 = c 2 and the "flow" is sonic. 



K<0 



K>0 



Figure 5. The Gaussian curvature on a torus: Doughnut surface or 
toroidal shell 



Therefore, the existence of an isometric immersion is equivalent to the existence of a 
weak solution of the balance equations of momentum in the fluid mechanic formulation, 
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which are nonlinear partial differential equations of mixed hyperbolic-elliptic type. Many 
usual surfaces have their Gauss curvature of changing sign, such as tori: Doughnut surfaces 
or toroidal shells in Fig. 

An appropriate approximate method to construct approximate solutions to the Gauss- 
Codazzi system has been designed, and a compensated compactness approach to establish 
the existence of a weak solution has been developed in [49j . The advantage of the compen- 
sated compactness approach is that it works for both the elliptic and hyperbolic phase. 
For more details about isometric embedding of Riemann manifolds in Euclidean spaces, 
we refer the reader to Han-Hong [83] and Yau [156] . 

For isometric embedding of a higher dimensional Riemannian manifold into the Eu- 
clidean space R with optimal dimension N (i.e., the Cartan- Janet dimension), the Gauss- 
Codazzi equations should be supplemented by the Ricci equations to describe the connec- 
tion form on the normal bundle. The Gauss-Codazzi-Ricci system even has no type in 
general, although it inherits important geometric features, including the beautiful Div- 
Curl structure which yields its weak continuity (see [50J). Also see Bryant-Griffiths- Yang 

na. 



5.2. Free-Boundary Techniques. To explain the techniques, we focus on the shock 
reflection-diffraction problem for potential flow which is widely used in aerodynamics. 

When a plane shock in the (t, x)-coordinates, x = (xi,X2) G R 2 , with the left state 
(/), V x $) = (pi, iti,0) and the right state (po,0,0),ui > 0, po < Pi, hits a symmetric 
wedge 

W := {x : |a?2| < x±taji9 w ,xi > 0} 

head on, it experiences a reflection-diffraction process, where p is the density and $ is the 
velocity potential of the fluid; see Fig. [TJ Then a self-similar reflected shock moves out- 
ward as the original shock moves forward in time. The complexity of reflection-diffraction 
configurations was first reported by Ernst Mach in 1878, and experimental, computational, 
and asymptotic analysis has shown that various different patterns may occur, including 
regular and Mach reflection (cf. [91 EH 11031 11161 1137} 1153] ). However, most of the fun- 
damental issues for shock reflection-diffraction have not been understood, including the 
global structure, stability, and transition of different patterns of shock reflection-diffraction 
configurations. Therefore, it becomes essential to establish the global existence and struc- 
tural stability of solutions in order to understand fully the shock reflection-diffraction 
phenomena. Furthermore, this problem is also fundamental in the mathematical theory of 
multidimensional conservation laws since their solutions are building blocks and asymp- 
totic attractors of general solutions to the two-dimensional Euler equations for compress- 
ible fluid flow (cf. Courant-Friedrichs [53], von Neumann [153j . and Glimm-Majda [81] : 
also see [HI EH1 EH 11031 IH6[ 1137] ) . As we will show below, the problem involves nonlin- 
ear partial differential equations of mixed elliptic-hyperbolic type, along with the other 
challenging difficulties such as free boundary problems and corner singularities, especially 
when a free boundary meets an elliptic degenerate curve. 

By scaling, the Euler equations for potential flow consist of the conservation law of mass 
and the Bernoulli law with the form: 



'd t p + V-(pV$) = 



a$ + i|v$i 2 + £-r = ^ 



,,-i oJ- (5-24) 



7— l 7— 1 ' 
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Shock wave 




Figure 6. Shock reflection-diffraction problem 

where 7 > 1 is the adiabatic exponent of the fluid. Then the reflection-diffraction problem 
can be formulated as the following mathematical problem. 

Initial-Boundary Value Problem. Seek a solution of the system with the initial 
condition at t = 0: 

(po,0) for \x%\ > Xitan9 w , x\ > 0, 

^(pi,uixi) for xi < 0, 

and the slip boundary condition along the wedge boundary dW : 

V$ • n\ dw = 0, (5.26) 

where n is the exterior unit normal to dW . 

Since the initial-boundary value problem is invariant under the self-similar scaling: 



(p,n 



t=0 



(5.25) 



(t,x)->(ai,ax), (p,$)->(p,-) 

a 

we seek self-similar solutions with the form: 



for a/0, 



for (£,77) = -. 



o(C 2 + V 2 ) is governed by the following 



Then the pseudo-potential function ip = ij) 
potential flow equation of second-order: 

div (p(|£V| 2 , tp)D<p) + 2p(\Dip\ 2 ,ip) = (5.27) 

with p(\Dtp\ 2 ,(p) := (pq~ - (7 + 1) ((p + i|D</?| 2 ))~, where D = (d^,d v ). The sonic speed 



is 



c(|£Vl ,y,Po 



7-11 



( P r i -(7-i)(^+-i^i 2 )) 



2n\V2 



(5.28) 
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Equation (|5.27p is a second-order nonlinear partial differential equations of mixed elliptic- 
hyperbolic type. It is strictly elliptic (i.e., pseudo-subsonic) if 

\D<p\ <c(\Dp\ 2 ,ip,p^ 1 ), (5.29) 

and strictly hyperbolic (i.e., pseudo-supersonic) if \Dcp\ > c(\D(p\ 2 , 92, /5q )• 

Shocks are discontinuities in the pseudo- velocity Dip. That is, if f2 + and il _ := $7 \ S7+ 
are two nonempty open subsets of f2 C M 2 and S := d£l + n fi is a C 1 curve where Dip has 
a jump, then ip G W^fi) n C 1 (il ± U S) n C 2 (n±) is a global weak solution of (l5\27D in 
in the sense of distributions if and only if cp is in Wi 'f°(fi) and satisfies equation (15.271) 
in Q 1 * 1 and the Rankine-Hugoniot condition on S: 

p(\D<p\ 2 ,<p)D<p -n\ s+ = p{\D<p\ 2 ,tp)D<p ■ n| g _. (5.30) 

The continuity of 92 is followed by the continuity of the tangential derivative of 92 across 
S. The discontinuity S of Dip is called a shock if ip further satisfies the physical entropy 
condition that the corresponding density function p(\Dip\ 2 , ip) increases across S in the 
pseudo-flow direction. 

The plane incident shock solution in the (£,x)-coordinates corresponds to a continuous 
weak solution ip in the self-similar coordinates (£, rj) with the following form: 

MC,r,) = -\(f + V 2 ) for £>£o, (5.31) 

M^v) = -\(f + r ] 2 )+u 1 (C-Co) for £< £o, (5.32) 

respectively, where £o > is the location of the incident shock, uniquely determined by 
G°o,/9i,7). Since the problem is symmetric with respect to the axis rj = 0, it suffices to 
consider the problem in the half-plane r/ > outside the half- wedge A := {£ < 0,7/ > 
0} U {rj > £ tan 0^, £ > 0}. Then the initial-boundary value problem can be formulated as 
the boundary value problem in A in the coordinates (£,77). 

Boundary Value Problem. Seek a solution ip of equation (I5.27P in the unbounded 
domain A with the slip boundary condition on the wedge boundary dA: 

Dip ■ n\ dA = (5.33) 

and the asymptotic boundary condition at infinity: 

JVo for £>&,«> £ tan 6 W , 2 2 

[<Pi for £ < to,V > 0, 
m i/te sense i/tat lim H92 — <^||c(A\b h (o)) = 0- 

Since <p\ does not satisfy the slip boundary condition, the solution must differ from 
<p\ in {£ < £0} H A, and thus a shock diffraction by the wedge occurs. A local existence 
theory of regular shock reflection-diffraction near the reflection point Pq can be established 
by following the von Neumann detachment criterion [153] with the structure of solution 
as in Fig. when the wedge angle is large and close to ^, in which the vertical line is 
the incident shock S = {£ = £0} that hits the wedge at the point Pq = (£o>£otan# w ), 
and state (0) and state (1) ahead of and behind S are given by 920 and 921, respectively. 
The solutions 92 and 921 differ only in the domain P0-P1-P2-P3 because of shock reflection 
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by the wedge boundary at Pq and diffraction by the wedge vertex P3, where the curve 
PqP\P2 is the reflected shock with the straight segment PqP\. State (2) behind PqP\ can 
be computed explicitly with the form: 

<£2(£,t?) = -tt(£ 2 + 7 ? 2 ) + ""2(£-£o) + (v ~ £0 t&n 6 w )u 2 t&n 6 W , (5.35) 



which satisfies Dip ■ n = 
PqPi are determined by 



on dA n {£ > 0}; the constant velocity «2 and the slope of 
&w, Po, Piil) from the two algebraic equations expressing (|5.3U|) 



and continuous matching of states (1) and (2) across PqP\, whose existence is exactly 
guaranteed by the condition on {6 W , po, Pi,l) which is necessary for the existence of a 
global regular shock reflection-diffraction configuration. 



Incident 
Shock 



(1) 
Sonic Circle 




RcficctC' 
Shock 



Figure 7. Regular reflection-diffraction configuration 



A rigorous mathematical approach has been developed in Chen-Feldman [38] to extend 
the local theory to a global theory for solutions of regular shock reflection-diffraction, 
which converge to the unique solution of the normal shock reflection when 9 W tends to 
5. The solution ip is pseudo-subsonic within the sonic circle for state (2) with center 
(«2j u 2 tan8 w ) and radius C2 > (the sonic speed) and is pseudo-supersonic outside this 
circle containing the arc -P1-P4 in Fig. [7J so that <p>2 is the unique solution in the domain 
P0P1P4. In the domain Q,, the solution is pseudo-subsonic, smooth, and C 1,1 -matching 
with state (2) across -P1-P4 and satisfies tp v = on P2P3; the transonic shock curve P1P2 
matches up to second-order with P0P1 continuously and is orthogonal to the £-axis at the 
point P2 so that the standard reflection about the £-axis yields a global solution in the 
whole plane. 

One of the main difficulties for the global existence is that the ellipticity condition for 
(J5.27P is hard to control. The second difficulty is that the ellipticity degenerates at the 
sonic circle -P1-P4. The third difficulty is that, on P1P4, the solution in Q has to be matched 
with ip2 at least in C 1 , i.e., the two conditions on the fixed boundary -P1-P4: the Dirichlet 
and conormal conditions, which are generically overdetermined for an elliptic equation 
since the conditions on the other parts of boundary are prescribed. Thus it is required 
to prove that, if (p satisfies (15.2TH in $7, the Dirichlet continuity condition on the sonic 
circle, and the appropriate conditions on the other parts of <9f2 derived from the boundary 
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value problem, then the normal derivative Dp ■ n automatically matches with Dp2 • n 
along -P1-P4. Indeed, equation (|5.27p . written in terms of the function w = p — p>2 in the 
(x, y)-coordinates defined near P1P4 such that -P1-P4 becomes a segment on {x = 0}, has 
the form: 

[2x — (7 + l)d x w)d xx w H — 2®vv w ~ ® xW = ® in 2; > and near x = 0, (5.36) 

C 2 

plus the "small" terms that are controlled by \ — 9 w in appropriate norms. Equation (|5.36p 
is strictly elliptic if d x w < — pj and become degenerate elliptic when d x w = — %■. Hence, 
it is required to obtain the C 1 ' 1 -estimates near -P1-P4 to ensure \d x w\ < -^ which in 
turn implies both the ellipticity of the equation in O and the match of normal derivatives 
Dip ■ n = D<p>2 ■ n along -P1-P4. Taking into account the "small" terms to be added to 
equation (|5.36j) . it is needed to make the stronger estimate \d x w\ < gT^xTj an d assume 
that ^ — 6 W is suitably small to control these additional terms. Another issue is the 
non- variational structure and nonlinearity of this problem which makes it hard to apply 
directly the approaches of Caffarelli [20J and Alt-Caffarelli- Friedman Od]. Moreover, the 
elliptic degeneracy and geometry of the problem makes it difficult to apply the hodograph 
transform approach in Kinderlehrer-Nirenberg [95] and Chen-Feldman [37] to fix the free 
boundary. 

For these reasons, one of the new ingredients in the approach is to develop further the 
iteration scheme in [36] to a partially modified equation. Equation (J5.36P is modified in Q 
by a proper Shiffmanization (i.e. cutoff) that depends on the distance to the sonic circle, 
so that the original and modified equations coincide for (p satisfying \d x w\ < „, jf , and 
the modified equation Mp = is elliptic in 0, with elliptic degeneracy on -P1-P4. Then a 
free boundary problem is solved for this modified equation: The free boundary is the curve 
P1P2, and the free boundary conditions on P1P2 are (p = <p\ and the Rankine-Hugoniot 
condition (|5.30[) . Moreover, the precise gradient estimate: \d x w\ < „/ §^ is made to 
ensure that <p satisfies the original equation (15.270 . 

This global theory for large-angle wedges has been extended in Chen-Feldman [39] up 
to the sonic angle 9 S < C , i.e. state (2) is sonic when 9 W = S , such that, as long 
as 9 W S (9 S ,^}, the global regular shock reflection-diffraction configuration exists, which 
solves the von Neumann's sonic conjecture (1943) [153j when u\ < c\. Furthermore, the 
optimal regularity of regular reflection-diffraction solutions near the pseudo-sonic circle 
has been shown in Bae-Chen-Feldman [6] to be C ' as established in [381139]. 

Some important efforts were made mathematically for the global reflection problem via 
simplified models, including the unsteady TSD equation and other nonlinear models (cf. 
|94|). Furthermore, in order to deal with the reflection problem, some asymptotic methods 
have been also developed including the work by Lighthill, Keller, Morawetz, and others. 
The physicality of weak Prandtl-Meyer reflection for supersonic potential flow around a 
ramp has been also analyzed in |66j. Another recent effort has been made on various 
important physical, mixed elliptic-hyperbolic problems in steady potential flow for which 
great progress has been made (cf . [281 E3 ESS 1117] and the references cited therein) . 

The self-similar solutions for the full Euler equations are required for the Mach reflection- 
diffraction configurations and general two-dimensional Riemann problems. For this case, 
the self-similar solutions are governed by a system that consists of two transport equations 
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and two nonlinear equations of mixed hyperbolic- elliptic type. One of the important fea- 
tures in the reflection-diffraction configurations is that the Euler equations for potential 
flow and the full Euler equations coincide in some important regions of the solutions. 

6. Singular Limits to Nonlinear Degenerate Hyperbolic Equations 

One of the important singular limit problems is the vanishing viscosity limit problem for 
the Navier-Stokes equations to the Euler equations for compressible barotropic fluids. The 
Navier-Stokes equations for a compressible viscous, barotropic fluid in Eulerian coordinates 
in M^_ take the following form: 

(d t p + d x (pv) = 0, 

\d t (pv) + d x (pv 2 + p) = ed xx v, 

with the initial conditions: 

p(0,x) = po(x), v(0,x)=v o (x) (6.2) 

such that lim. x ^.± 00 (pQ(x),vo(x)) = (/3 ± ,t> ± ), where p denotes the density, v represents 
the velocity of the fluid when p > 0, p is the pressure, m = pv is the momentum, and 
(/3 ± ,u ± ) are constant states with p^ > 0. The physical viscosity coefficient e is restricted 
toe £ (0, so] for some fixed e$ > 0. For a polytropic perfect gases, the pressure-density 
relation is determined by (|3.14p - (|3.15p with adiabatic exponent 7 > 1. 

Formally, when e — > 0, the Navier-Stokes equations become the isentropic Euler equa- 
tions (|3.13p . which is strictly hyperbolic when p > 0. However, near the vacuum p = 0, the 
two characteristic speeds of (|3. 13[) may coincide and the system be degenerate hyperbolic. 

The vanishing artificial/numerical viscosity limit to the isentropic Euler equations with 
general L°° initial data has been studied by DiPerna [61J, Chen [29, 31J, Ding [58j, 
Ding-Chen-Luo |59j . Lions- Perthame-Souganidis [106J, and Lions-Perthame-Tadmor [107J 
via the methods of compensated compactness. Also see DiPerna [62], Morawetz [118j, 
Perthame-Tzavaras [129] . and Serre |138j for the vanishing artificial/numerical viscosity 
limit to general 2x2 strictly hyperbolic systems of conservation laws. The vanishing 
artificial viscosity limit to general strictly hyperbolic systems of conservation laws with 
general small BV initial data was first established by Bianchini-Bressan |15j via direct BV 
estimates with small oscillation. Also see LeFloch-Westdickenberg [104| for the existence 
of finite-energy solutions to the isentropic Euler equations with finite-energy initial data 
for the case 1 < 7 < 5/3. 

The idea of regarding inviscid gases as viscous gases with vanishing real physical vis- 
cosity can date back the seminal paper by Stokes [145J and the important contribution 
of Rankine [132], Hugoniot [85], Rayleigh [155] . and Taylor [HH] (cf. Dafermos [53]). 
However, the first rigorous convergence analysis of vanishing physical viscosity from the 
Navier-Stokes equations I16.1J) to the isentropic Euler equations was made by Gilbarg [77] 
in 1951, when he established the mathematical existence and vanishing viscous limit of 
the Navier-Stokes shock layers. For the convergence analysis confined in the framework of 
piecewise smooth solutions; see Hoff-Liu [87], Gues-Metivier-Williams-Zumbrun [82] . and 
the references cited therein. The convergence of vanishing physical viscosity with general 
initial data was first studied by Serre-Shearer [139] for a 2 x 2 system in nonlinear elasticity 
with severe growth conditions on the nonlinear function in the system. 
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In Chen-Perepelitsa [15], we have first developed new uniform estimates with respect 
to the real physical viscosity coefficient for the solutions of the Navier-Stokes equations 
with the finite-energy initial data and established the if -1 -compactness of weak entropy 
dissipation measures of the solutions of the Navier-Stokes equations for any weak entropy 
pairs generated by compactly supported C 2 test functions. With these, the existence of 
measure-valued solutions with possibly unbounded support has been established, which 
are confined by the Tartar-Murat commutator relation with respect to two pairs of weak 
entropy kernels. Then we have established the reduction of measure- valued solutions with 
unbounded support for the case 7 > 3 and, as corollary we have obtained the existence 
of global finite-energy entropy solutions of the Euler equations with general initial data 
for 7 > 3. We have further simplified the reduction proof of measure- valued solutions 
with unbounded support for the case 1 < 7 < 5/3 in LeFloch-Westdickenberg [104j and 
extended to the whole interval 1 < 7 < 3 . With all of these, we have established the 
first convergence result for the vanishing physical viscosity limit of solutions of the Navier- 
Stokes equations to a finite-energy entropy solution of the isentropic Euler equations with 
general finite-energy initial data. We remark that the existence of finite-energy solutions 
to the isentropic Euler equations for the case 7 > 5/3 have been established, which is in 
addition to the result in [104J for 1 < 7 < 5/3. 

More precisely, consider the Cauchy problem f)6.ip — f|6.2j) for the Navier-Stokes equations 
in R^_. Hoff's theorem in [86] (also see Kanel [90J for the case of the same end-states) 
indicates that, when the initial functions (po(x),vo(x)) are smooth with the lower bounded 
density po(x) > eg > for x £ R and 

lim (p (x),v (x)) = (p ± ,v ± ), 

x— >±oo 

then there exists a unique smooth solution (p £ (t, x),v £ (t, x)), globally in time, with p £ (t, x) > 
c £ (t) for some c £ (t) > for t > and lim x _ ! ,± 00 (p £ (t, x),v £ (t, x)) = (p^jV^). 

Let (p(x),v(x)) be a pair of smooth monotone functions satisfying (p(x),v(x)) = (p^^^) 
when ±x > Lq for some large L$ > 0. The total mechanical energy for (J6.1D in R with 
respect to the pair (p, v) is 

E[p,v](t) := / (r]*(p,rn)—r]*(p,m) — Vri*(p 1 fh)-(p — p,rn — fh))dx>0, (6.3) 

where fh = pv. 

Combining the uniform estimates and the compactness of weak entropy dissipation 
measures in H7~ with the compensated compactness argument and the reduction of the 
measure- valued solution f i)X , we conclude 

Theorem 6.1 (Chen-Perepelitsa [45]). Let the initial functions (pq,Vq) be smooth and 
satisfy the following conditions: There exist Eq,Ei,Mq > 0, independent of e, and Cq > 
such that 

(i) pl(x) > c £ > 0, j P £ (x)\v £ (x) - v(x)\dx < M < 00; 
(ii) The total mechanical energy with respect to (p, u) is finite: 

/ (2 Po^) K 0*0 " v(x)\ 2 + h(p e (x),p(x))jdx <E <oo, 

where h(p, p) = pe(p) - pe(p) - (e(p) + pe'(p)) (p - p); 
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(iii) e 2 r JfU*)!! d < E 

(iv) (pq(x), Pq(x)vq(x)) — > (po(x) , po(x)vo(x)) in the sense of distributions as e — > 0, 
with po{x) > a.e. 

Let (p £ ,m £ ),m £ = p £ v £ , be the solution of the Cauchy problem (I6.ip - (|6,2p for the Navier- 
Stokes equations with initial data (pq(x) , pf > (x)vQ(x)) for each fixed e > 0. Then, when 
e — > 0, there exists a subsequence of (p £ ,m £ ) that converges almost everywhere to a finite- 
energy entropy solution (p, m) of the Cauchy problem (J3.13P and (16.2f) with initial data 
(po(x), pq(x)vq(x)) for the isentropic Euler equations with 7 > 1. 
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